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UNDER THE T.TLE 


HIGHER MATHEMATICS. 


First Edition, September, 1896. 
Second Edition, January, 1898. 
Third Edition, August, 1900. | 
Fourth Edition, January, 1906, 


BOLT ORS RUREGACE, 


Tur volume called Higher Mathematics, the first edition 
of which was published in 1896, contained eleven chapters by 
eleven authors, each chapter being independent of the others, 
but all supposing the reader to have at least a mathematical 
training equivalent to that given in classical and engineering 
colleges. The publication of that volume is now discontinued 
and the chapters are issued in separate form. In these reissues 
it will generally be found that the monographs are enlarged 
by additional articles or appendices which either amplify the 
former presentation or record recent advances. This plan of 
publication has been arranged in order to meet the demand of 
teachers and the convenience of classes, but it is also thought 
that it may prove advantageous to readers in special lines of 
mathematical literature. 

It is the intention of the publishers and editors to add other 
monographs to the series from time to time, if the call for the 
same seems to warrant it. Among the topics which are uncer 


- consideration are those of elliptic functions, the theory of num- 


bers, the group theory, the calculus of variations, and non- 
Euclidean geometry; possibly also monographs on branches of 
astronomy, mechanics, and mathematical physics may be included. 
It is the hope of the editors that this form of publication may 
tend to promote mathematical study and research over a wider 
field than that which the former volume has occupied. 


December, 1905. * 
iii 


AUTHOR’S PREFAGE: 


In republishing this short treatise in book form the author 
solicits criticism but offers no apology. The type of the book 
he has sought to imitate is that shown in the “ mathematical 
tracts” of the late Sir George B. Airy. The brevity and the con- 
crete illustrations of these “ tracts”? have served very effectively 
in introducing students to a number of the more difficult fields 
of applied mathematics; and it is hoped that this treatise will 
serve a similar end. 

The theory of probability and the theory of errors now con- 
stitute a formidable body of knowledge of great mathematical 
interest and of great practical importance. Though developed 
largely through applications to the more precise sciences of as- 
tronomy, geodesy, and physics, their range of applicability extends 
to all of the sciences; and they are plainly destined to play an 
increasingly important réle in the development and in the appli- 
cations of the sciences of the future. Hence their study is not 
only a commendable element in a liberal education, but some 
knowledge of them is essential to a correct understanding of 
daily events. 

No special novelty of presentation is claimed for this work; 
but the reader may find it advantageous to know that a definite 
plan has been followed. This plan consists in presenting each 
principle, first, by means of a simple, concrete example; passing, 
secondly, to a general statement by means of a formula; and, 
thirdly, illustrating applications of the formula by concrete 
examples. Great pains have been taken also to secure clear and 
correct statements of fundamental facts. If these latter are 
duly understood, the student needs little additional aid; if 
they are not duly understood, no amount of aid will forward him. 

The passage from the elementary concrete to the advanced 
abstract may appear to be abrupt to the reader in some cases. 
It is hoped, however, that any large gaps may be easily bridged 
and that any serious difficulties may be easily overcome by means 
of the references given to the literature of the subject. In any 
event the student will find that in this, as in all of the more ardu- 
ous sciences, his greatest pleasure and his highest discipline will 
~ come from bridging such gaps and from surmounting such diffi- 

culties. 

WasHINGTON, D. C., December, 1905. 
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PROBABILITY AND THEORY OF ERRORS. 


ART. 1. INTRODUCTION. 


IT is a curious circumstance that a science so profoundly 
mathematical as the theory of probability should have origi- 
nated in the games of chance which occupy the thoughtless and 
the profligate.** That such is the case is sufficiently attested 
by the fact that much of the terminology of the science and 
many of its familiar illustrations are drawn directly from the 
vocabulary and the paraphernalia of the gambler and the trick- 
ster. It is somewhat surprising, also, considering the antiquity 
of games of chance, that formal reasoning on the simpler 
questions in probability did not begin before the time of Pascal 
and Fermat. Pascal was led to consider the subject during the 
year 1654 through a problem proposed to him by the Chevalier 
de Méré, a reputed gambler.t The problem in question is 
known as the problem of points and may be stated as follows: 
two players need each a given number of points to win ata 
certain stage of their game; ifthey stop at this stage, how should 
the stakes be divided? Pascal corresponded with his friend 
Fermat on this question; and it appears that the letters which 
passed between them contained the earliest distinct formulation 
of principles falling within the theory of probability. These 


* The historical facts referred to in this article are drawn mostly from Tod- 
hunter’s History of the Mathematical Theory of Probability from the time of 
Pascal to that of Laplace (Cambridge and London, 1865). 

+‘ Un probléme relatif aux jeux de hasard, proposé aun austére janséniste 
par un homme du monde, a été l’origine du calcul des probabilités.” Poisson, 
Recherches sur la Probabilité des Jugements (Paris, 1837). 
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acute thinkers, however, accomplished little more than a correct 
start in the science. Each seemed to rest content at the time 
with the approbation of the other. Pascal soon renounced 
such mundane studies altogether; Fermat had only the scant 
leisure of a life busy with affairs to devote to mathematics; 
and both died soon after the epoch in question,—Pascal in 
1662, and Fermat in 1665. 

A subject which had attracted the attention of such dis- 
tinguished mathematicians could not fail to excite the interest 
of their contemporaries and successors. Amongst the former 
Huygens is the most noted. He has the honor of publishing 
the first treatise* on the subject. It contains only fourteen 
propositions and is devoted entirely to games of chance, but it 
gave the best account of the theory down to the beginning of 
the eighteenth century, when it was superseded by the more elab- 
orate works of James Bernoulli,t Montmort,t and De Moivre.§ 
Through the labors of the latter authors the mathematical 
theory of probability was greatly extended. They attacked, 
quite successfully in the main, the most difficult problems; 
and great credit is due them for the energy and ability dis- 
played in developing a science which seemed at the time to 
have no higher aim than intellectual diversion. Their names, 
undoubtedly, with one exception, that of Laplace, are the most 
important in the history of probability. : 

Since the beginning of the eighteenth century almost every 
mathematician of note has been a contributor to or an expos- 
itor of the theory of probability. Nicolas, Daniel, and John 
Bernoulli, Simpson, Euler, d’Alembert, Bayes, Lagrange, Lam- 
bert, Condorcet, and Laplace are the principal names which 
figure in the history of the subject during the hundred years 

* De Ratiociniis in Ludo Alee, 1657. 
+ Ars Conjectandi, 1713. 


t Essai d’Analyse sur les Jeux de Hazards, 1708. 
§ The Doctrine of Chances, 1718. 


| Todhunter says of Montmort, for example, “In 1708 he published his 


work on Chances, where with the courage of Columbus he revealed a new world 
to Mathematicians,” 
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ending with the first quarter of the nineteenth century. Of 
contributions from this brilliant array of mathematical talent, 
the Théorie Analytique des Probabilités of Laplace is by far 
the most profound and comprehensive. It is, like his Mé- 
canique Céleste in dynamical astronomy, still the most elabo- 
rate treatise on the subject. An idea of the grand scale of the 
work in its present form* may be gained by the facts that the 
non-mathematical introductiont+ covers about one hundred and 
fifty quarto pages; and that, in spite of the extraordinary 
brevity of mathematical language, the pure theory and its ac- 
cessories and applications require about six hundred and fifty 
pages. 

From the epoch of Laplace down to the present time the 
extensions of the science have been most noteworthy in the 
fields of practical applications, as in the adjustment of obser- 
vations, and in problems of insurance, statistics, etc. Amongst 
the most important of the pioneers in these fields should 
be mentioned Poisson, Gauss, Bessel, and De Morgan. Nu- 
merous authors, also, have done much to simplify one or an- 
other branch of the subject and thus bring it within the range 
of elementary presentation. The fundamental principles of 
the theory are, indeed, now accessible in the best text-books 
onalgebra: and there are many excellent treatises on the pure 
theory and its various applications. 

Of all the applications of the doctrine of probability none 
is of greater utility than the theory of errors. In astronomy, 
geodesy, physics, and chemistry, as in every science which at- 
tains precision in measuring, weighing, and computing, a 
knowledge of the theory of errors is indispensable. By the aid 
of this theory the exact sciences have made great progress dur- 


*The form of the third edition published in 1820, and of Vol. VII of the 
complete works of Laplace recently republished under the auspices of the 
Académie des Sciences by Gauthier-Villars. This Vol. VII bears the date 1886. 

+ ‘* Cette Introduction,” writes Laplace, ‘‘est le developpement d’une Legon 
sur les Probabilités, que je donnai en 1795, aux Ecoles Normales, oti je fus ap- 
pelé comme professeur de Mathématiques avec Lagrange, par un décret de la 

_ Convention nationale.” 
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ing the nineteenth century, not only in the actual determination 
of the constants of nature, but also in the fixation of clear 
ideas as to the possibilities of future conquests in the same di- 
rection: Nothing, for example, is more satisfactory and in- 
structive in the history of science than the success with which 
the unique method of least squares has been applied to the 
problems presented by the earth and the other members of the 
solar system. So great, in fact, are the practical value and 
theoretical importance of the method of least squares, that it is 
frequently mistaken for the whole theory of errors, and is 
sometimes regarded as embodying the major part of the doc- 
trine of probability itself. 

As may be inferred from this brief sketch, the theory of 
probability and its more important applications now constitute 
an extensive body of mathematical principles. and precepts. 
Obviously, therefore, it will be impossible within the limits of 
a single condensed monograph to do more than give an out- 
line of the salient features of the subject. It is hoped, how- 
ever, in accordance with the general plan of the volume, that 
such outline will prove suggestive and helpful to those who 
may come to the science for the first time, and also to those 
who, while somewhat familiar with the difficulties to be over- 
tome, have not acquired a working knowledge of the subject. 
Effort has been made especially to clear up the difficulties of 
the theory of errors by presenting a somewhat broader view of 
the elements of the subject than is found in the standard 
treatises, which confine attention almost exclusively to the 
method of least squares. This chapter stops short of that 
method, and seeks to supply those phases of the theory which 
are either notably lacking or notably erroneous in works 
hitherto published. It is believed, also, that the elements here 
presented are essential to an adequate understanding of the 
-well-worked domain of least squares.* 


x . . , : i 
The author has given a brief but comprehensive statement of the method 


_ of least squares in the volume of Geographical Tables published by the Smith- 
sonion Institution, 1894. 
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ART. 2. PERMUTATIONS. 

The formulas and results of the theory of permutations 
and combinations are often needed for the statement and so- 
lution of problems in probabilities. This theory is now to be 
found in most works on algebra, and it will therefore suffice 
here to state the principal formulas and illustrate their mean- 
ing by a few numerical examples. 

The number of permutations of z things taken 7 in a group 
is expressed by the formula 

(x), = nw —1)'n — 2)... (#—7r+1). (1) 

Thus, to‘illustrate, the number of ways the four letters a, b, 
c,d can be arranged in groups of two is 4. 3 = 12, and the groups 
are ; 

Gb. bas eimta, ad), da, be, cb, Od, gb, cay: ges 

Similarly, the formula gives for 

ie and y= 2, (3) =r 3}..2 =.0, 

ha a PER (7) 07.22 s=I2 10, 

m.=10 -.n=6, (10),= 10.9.8.7.6.'5 = 151200. 


The results which follow from equation (1) when and r 
do not exceed 10 each are embodied in the following table: 


VALUES OF PERM UTATIONS. 


“1814400 | 362880 40320 
3628800 | 362880 
3628800 


10 9 8 7 6 5;) 44-3421 1 

I 10 ax) 8 f 4 6 LA age el each allt Cabo fe 
2 go 72 £6 7a 1°) 20 | 12 Onie2 

3 720 BOs eens oO 210 120 60 | 24 6 

‘4 5040 3024 1680 |- 840 360 | 120 | 24 

5. 302.40, 15120 6720 | 2520 FZON E20. py 

ii 6 ' “151200 60480 20160 5040 720 : 

7. 604800 | 181440 40320 5040 

3 

9 
te) 


Lal 


Sp 9864100 | 986409 | 109600 13699 | 1956 | 325 | 64] 15 | 4] I 


The use of this table is obvious. Thus, the number of per- 
mutations of eight things in. groups of five each is found in the 
fifth line of the column headed with the number 8. It will be 
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noticed that the last two numbers in each column (excepting 
that headed with 1) are the same. This accords with the for- 
mula, which gives for the number of permutations of 2 things 
in groups of the same value as for ” things in groups of 
(7 —1). It will also be remarked that the last number in each 
column of the table is the factorial, ~!, of the number x at the 
head of the column, For example, in the column under 7, the 
last nufiber is 5049 = 14223. 45) 0n7 — 74 

The total number of permutations of x things taken singly, 
in groups of two, three, etc., is found by summing the numbers 
given by equation (1) for all values of 7 from 1 to.” Calling 
this total or sum S,, it will be given by 

Ss = 2(72),« (2) 

To illustrate, suppose z = 3, and, to fix the ideas, let the 
three things be the three digits 1, 2, 3. Then from the above 
table it is seen that S, = 3 -+6+6= 15; or, that the number 
of numbers (all different) which can be formed from those dig- 
its is fifteen. These numbers are 1, 2, 35 12,13, 21,2354 teoo- 
128, 132, 213, 231, 312, 32% 

The values of S, for z= 1, 2,,..10 are given under the 
corresponding columns of the above table. But when x is 
large the direct summation indicated by (2) is tedious, if not 
impracticable. Hence a more convenient formula is desirable. 
To get this, observe that (1) may be written 

a! 

(= Eo, (1) 
if 7 is restricted to integer values between 1 and (x — 1), both 
inclusive. This suffices to give all terms which appear in the 
right-hand member of (2), since the number of permutations 


for 7 = (” — 1) is the same as forr =~. Hence it appears 
that 
Ah me! n!\ 


S,=n2!+—+——4,,. 


I I.2 
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But as increases, the series by which 7! is here multiplied 
approximates rapidly towards the base of natural logarithms ; 
that is, towards 

e = 2.7182818 +, log ¢ = 0.4342945. 
Hence for large values of x 
5 = n!e, approximately.* (3) 
To get an idea of the degree of approximation of (3), sup- 


pose z=g. Then the computation runs thus (see values in 
the above table): 


log 
g! = 362880 5.5597630 
é 0.43429045 


g!e= 986410 —_§.9940575 
S, = 986409 by equation (2). 
The error in this case is thus seen to be only one unit, or 
about one-millionth of S,.+ 


Prob. r. Tabulate a list of the numbers of three figures each 
which can be formed from the first five digits 1,...5. Howmany 
numbers can be formed from the nine digits? 

Prob. 2. Is S, always an odd number for 2 odd? Observe 
values of S;, in the table above. 


ART. 3. COMBINATIONS. 


In permutations attention is given to the order of arrange- 
ment of the things considered. In combinations no regard is 
paid to the order of arrangement. Thus, the permutations of 
the letters a, 4, c, din groups of three are 


(abc) (abd) bac bad ach (acd) cab cad 
adb- ade dab dac bea (bcd) cha chad 
bda bde dba dbe cda cab ada aeb 


* See Art. 6 for a formula for computing #! when x is a large number. 

+ When large numbers are to be dealt with, equations (1)' and (3) are easily 
managed by logarithms, especially if a table of values of log (7!) is available. 
Such tables are given to six places in De Morgan’s treatise on Probability in 
the Encyclopedia Metropolitana, and to five places in Shortrede’s Tables 


(Vol. I, 1849). 


14 PROBABILITY AND THEORY OF ERRORS. 


But if the order of arrangement is ignored all of these are 
seen to be repetitions of the groups enclosed in parentheses, 
namely, (abc), (abd), (acd), (bcd). Hence in this case out of 
twenty-four permutations there are only four combinations. 


A general formula for computing the number of combina- 
tions of 2 things taken in groups of 7 things is easily derived. 
For the number of permutations of things in groups of 7 is 
by (1) of Art. 

ane = nn —'Ilnv—2)...(@—r+ 1); 
and since each group of 7 things gives 1.2.3...7= 7! per- 
mutations, the number of combinations must be the quotient 
of (z), by !. Denote this number by C(z),. Then the gen- 
eral formula is 
Gee n(z— 1)(#—2)...(2—7r-+1) 


7! 


(1) 


This formula gives, for example, in the case of the four let- 
ters a, 6, c,d taken in groups of three, as considered above, 


42 
C(4), = fee = 4. 
Multiply both numerator En denominator of the right-hand 
member of (1) by (z — 7)! The result is 
n!| 


C(n), = ri(z—nv (1) 


which shows that the number of combinations of z things in 
groups of y is the same as the number of combinations of z 
things in groups of (z — 7). . Thus, the number of combina- 
tions of the first ten letters a, 6,c...y in groups of three or 


seven is 
10! 
3! 317! = 120. 
The following table gives the values C(x), for all values of 
mand r from I to fo. 


The mode of using this table is evident. For example, the 


number of combinations of eight things in sets of five each is | 


found on the fifth line of the column headed 8 to be 50. 
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VALUES OF COMBINATIONS. 


ROD 9 8 FPG Sat a 3) 2% 
I 10 9 8 7 6 5 4|3 I 
2 45 36 28 2 Tats fo) Ol saior 
ei 120 84 56 a5) 20 0 tO AE 
4 210 | 126 70 Ch Wily 5 I 
5 252 | 126 56 21 6 I 
6 210 84 28 7 I 
7 120 36 8 I 
8 45 9 I 
9 10 I 
10 I 
Se ROSS ues She eos selel yal OS nes ballets wens les 


It will be observed that the numbers in any column show 
a maximum value when z is even and two equal maximum 
values when z is odd. That this should be so is easily seen 
from (1)’, which shows that C(7), will be a maximum for any 
value of 2 when 7!(z — 7)! isa minimum. For x even this is 
a minimum for 7 = 42; while for z odd it has equal minimum 
values for 7 =4(z — 1) andr=43(u-+1). Thus, 


: na} 
maximum of C(m), = 7— , for # even, 
) 


(2 
ze 
n' 


Tepe PU al 
Demin RO. 


(2) 


for z odd. 


The total number of combinations of z things taken singly, 
in groups of two, three, etc., is found by summing the numbers 
given by (1) for all values of 7 from I to both inclusive. 
Calling this total or sum Sis 

Sears Ce 
The same sum will also come from (1)’ by giving to 7 all values 
from 1 to (7 — 1), both inclusive, summing the results, and in- 
creasing their aggregate by unity. Thus by either process 


n(z — 1) 


san MR OD pat 


125;3 
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The second member of this equation is evidently equal to 
(1+1)*—1. Hence 

Sp 2C(aje sre ae (3) 
The values of S, for values of x and v from I to Io are given 
under the corresponding columns of the above table. 


Prob. 3. How many different squads of ten men each can be 
formed from a company of 100 men? 

Prob. 4 How many triangles are formed by six straight lines 
each of which intersects the other five ? 

Prob. 5. Examine this statement: “In dealing a pack of cards 
the number of hands, of thirteen cards each, which can be produced 
is 635 013 559 600. But in whist four hands are simultaneously held, 
and the number of distinct deals . . . would require twenty-eight 
figures to express it.” * 


Prob. 6. Assuming combination always possible, and disregarding 
the question of proportions, find how many different substances 
could be produced by combining the seventy-three chemical ele- 
ments. 


ART. 4. DIRECT PROBABILITIES. 


If it is known that one of two events must occur in any 
trial or instance, and that the first can occur in @ ways and the 
second in 6 ways, all of which ways are equally likely to hap- 
pen, then the probability that the first will happen is expressed 
mathematically by the fraction a/(a + 4), while the probability 
that the second will happen is 6/(a-+- 6). Such events are said 
to be mutually exclusive. Denote their probabilities by p and 
g respectively. Then there ee 


a 
hE ha ee Page (1) 
the last equation following from the first two and being the 
mathematical expression for the certainty that one of the two 
events must happen. 
Thus, to illustrate, in tossing a coin it must give “head ” or 
“tail”; a= b=1, and p= g=1/2. Again, if an utn contain 
a =5 white and 6 =8 black balls, the probability of drawing 


* Jevons, Principles of Science, New York, 1874, p. 217. 


ee 
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a white ball in one trial is f= 5/13 and that of drawing a 
black one g = 8/13. 

Similarly, if there are several mutually exclusive events 
which can occur in a, 6, ¢... ways respectively, their probabil- 
ies Dg. 7 a>. are giver by 

a b c 
ee eh ht a bac a,’ eee 
(2 


2teatrt...=. 


For example, if an urn contain a= 4 white, d= 5 black, 


and c = 6 red balls, the probabilities of drawing a white, black, 
and red ball at @ siecle trial are p= 4/15, g = $/15, and 
r = 6/15, respectively. 

Formulas (1) and (2) may be applied to a wide variety of 
cases, but it must suffice here to give only a few such. As a 
first illustration, consider the probability of drawing at random 
a number of three figures from the entire list of numbers which 
can be formed from the first seven digits. A glance at the 
table of Art. 1 shows that the symbols of formula (1) have in 
this case the values a= 210, and a+6= 13699. Hence 
b = 13489, and fp = 210/13699; that is, the probability in ques- 
tion is about 1/65. 

Secondly, what is the probability of holding in a hand of 
whist all the cards of one suit? Formula (1) of Art. 3 shows 
that the number of different hands of thirteen cards each which 
may be formed from a pack of fifty-two cards is 

B2UST BO... 40 

ae ye ere Be 253.559 000, 
and the probability required is the reciprocal of this number. 
The probability against this event is, therefore, very nearly 
unity. 

Thirdly, consider the probabilities presented by the case of 
an urn containing 4 white, 5 black, and 6 red balls, from which 
at a single trial three balls are to be drawn. Evidently the 
triad of balls drawn may be all white, all black, all red, partly 
white and black, partly white and red, partly black and red, or 
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one each of the white, black, and red. There are thus seven 
different probabilities to be taken into account. The theory 
of combinations shows (see equation (1), Art. 3) that the total 


number of 


White triads = = 4=a 
T s2.053 
Black triads =o = 10=5 
Red triads = z 3 ae ane Oh=n Cc 
White and black triads eh Pot ahy —( 4+10)= jo=d 
= == 4 7 
White and red triads = 2 Be —( 4+20)== 96 =e 
ee Pig 10.69 
Black and red triads eS —¢ = (10-+-20) =135 =the 
White, black, and red triads = 4.5.6 = 120 ee 


Sum = 455 


The total number of these triads is 455, and is, as it should 
be, the number of combinations in groups of three each of the | 
whole number of balls. Hence formulas (2) give the seven 
different probabilities which follow, using the initial letters 


w, 6, r to indicate the colors represented in a triad: 


For a triad zw P=) AAs: 
ae oe 12) q= 10/455, 
eee eg eo =O Fae e 
“eo 6 ww or whb So) OLAS 5, 
(co 6 wir or wre je GO) AGG. 
Be OOP Ow brn “2135/4585; 
CEO. tube v = 120/455. 


Prob. 7. When three dice are thrown together, what is the prob- 


ability that the throw will be greater than g? 


Prob. 8. Write down a literal formula for the probabilities of the 
several possible triads considered in the above question of the balls, 
supposing the numbers of white, black, and red balls to be Z, m, n, 


respectively. 
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ArT. 5. PROBABILITY OF CONCURRENT EVENTS. 


If the probabilities of two independent events are f, and 
pf. respectively, the probability of their concurrence in any 
single instance is ~,f, Thus, suppose there are two urns 
U, and U,, the first of which contains @, white. and 4, black 
balls, and the second a, white and 4, black balls. Then the 
probability of drawing a white ball from UV, is ~,= a,/(a, + 4,), 
while that of drawing a white ball from U, is ~, = a,/(a, + 4). 
The total number of different pairs of balls which can be formed 
from the entire number of balls is (a, + ,)(@,-+4,). Of these 
pairs a,a, are favorable to the concurrence of white in simul- 
taneous or successive drawings from the two urns. Hence the 
- probability of a concurrence of 

white with white = a,a,/(a, + 6,)(@, + 4), 

white with black = (@,4,+ a,0,)/(a, + 4,)(@, + 4), 

black with black = 6,4,/(a, + 6,)(@, + 4), 
and the sum of these is unity, as required by equations (2) of 
Article 4. 

In general, if 2,, 2,, ,.- - denote the probabilities of several 
independent events, and P denote the probability of their 
concurrence, 

Rise PIPE es (1) 

To illustrate this formula, suppose there is required the 
probability of getting three aces with three dice thrown simul. 
taneously. In this case ~, =f, =p, = 1/6 and 
Ph /6), = 1/210. 

Similarly, if two dice are thrown simultaneously the proba- 
bility that the sum of the numbers shown will be 11 is 2/36; 
and the probability that this sum II will appear in two succes- 
sive throws of the same pair of dice is 4/36.36. 

The probability that the alternatives of a series of events 
will concur is evidently given by 

Q= 999s -++ = (0 Pit — p(t 7a) chen © (2) 

Thus, in the case of the three dice mentioned above, the 

probability that each will show something other than an ace is 
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9, = 9, =% = 5/6, and the probability that they will concur in 
thissis Q)==1125,/210. 

Many cases of interest occur for the application of (1) and 
(2). One of the most important of these is furnished by suc< 
cessive trials of the same event. Consider, for example, what 
may happen in z trials of an event for which the probability 
is p and against which the probability is g. ‘The probability 
that the event will occur every time is f”. The probability that 
the event will occur (z — 1) times in succession and then fail is 
p’-'q. But if the order of occurrence is disregarded this last 
combination may arrive in # different ways; so that the prob. 
ability that the event will occur (z — 1) times and fail once is 
np"~'g. Similarly, the probability that the event will happer 
(2 — 2) times and fail twice is $u(z — 1) p'~*g'; etc. That ts, 
the probabilities of the several possible occurrences are given by 
the corresponding terms in the development of (f+ @q)”. 


By the same reasoning used to get equations (2) of Art. 
3 it may be shown that the maximum term in the expansion 
of (p+ ¢)” is that in which the exponent m, say, of g is 
the whole number lying between (%~-+ 1)g — 1 and (z+ 1)g¢. 
In other words, the most probable result in z trials is the 
occurrence of the event (z—m) times and its failure m 
times. When ~ is large this means that the most probable of 
all possible results is that in which the event occurs 2 — ng 
= n(1 — g) = up times and fails zgtimes. Thus, if the event 
be that of throwing an ace with a single die the most probable 
of the possible results in 600 throws is that of 100 aces and 
500 failures. 

Since g” is the probability that the event will fail every time 
in # trials, the probability that it will occur at least once in z 
trials is 1 —g”. Calling this probability 7,* 


Ae Sak Bg 8a) (3) 
If v in this equation be replaced by 1/2, the corresponding 
value of z is the number of trials essential to render the 


* See Poisson’s Probabilité des Jugements, pp. 40, 41. 
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chances even that the event whose probability is f will occur 
at least once. Thus, in this case, the value of z is given by 
log 2 

log (1—p) 

This shows, for example, if the event, be the throwing of double 
sixes with two dice, for which » = 1/36, that the chances are 
even (7 = 1/2) that in 25 throws (z = 24.614 by the formula) 
double sixes will appear at least once. 


xu — 


Equation (3) shows that however small # may be, so long as 
it is finite, 7 may be taken so large as to make 7 approach in- 
definitely near to unity; that is, 2 may be so large as to render 
it practically certain that the event will occur at least once. 


When ~ is large 


n(n” — 1) » a(n —I1)\("— 2), 


C7) Se ae a p+... 


Seeds ve ~ a 3 ap 


= e-"* approximately. 


Thus an approximate value of 7 is 


roi eo log ¢ = 0.4342495. (4) 
‘This formula gives, for example, for the probability of drawing 
the ace of spades from a pack of fifty-two cards at least once in 
104 trials y= 1 — e~* = 0.865, while the exact formula (3) 
gives 0.867. 

Similarly, the probability of the occurrence of the event at 
least ¢ times in x trials will be given by the sum of the terms 
of (p+ @)” from fp” up to that in f’g”~‘ inclusive. This proba- 
bility must be carefully distinguished from the probability that 
the event will occur ¢ times only in the z trials, the latter being 


expressed by the single term in pg". 


_ Prob. 9. Compare the probability of holding exactly four aces in 
five hands of whist with the probability cf nolding at least four aces 
in the same number of hands. 

Prob. 10. What is the probability of an event if the chances are 
even that it occurs at least once in a million trials? See equation (4). 
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ArT. 6. BERNOULLI’S THEOREM. 


Denote the exponents of f and g in the maximum term of 
(p+ g)” by wand m respectively, and denote this term by 7. 
Then 


pe eee CR 2) eee a! 


i mn =_— 
m | fee im! 


pig”. (1) 

As shown in Art. 5, # in this formula is the greatest whole 
number in (z+ 1)p, and m the greatest whole number in 
(x -+ 1)g; so that when » is large, wand m are sensibly equal 
to zp and ng respectively. 

The direct calculation of 7 by (1) is impracticable when z 
is large. To overcome this difficulty the following expression 
is used ;* 


ade aii I 
ni=n"e-” Vann t Sit ae alee $ ). (2) 


12” 


log e = 0.4342495, log 27 =0.7981799. 


This expression approaches "e-" V¥22” as a limit with the 
increase of z, and in this approximate form is known as Stir- 
ling’s theorem. Although a rude approximation to #! for 
small values of z this theorem suffices in nearly all cases 
wherein such probabilities as Z are desired. Making use of 
the theorem in (1) it becomes 


Se Eras (3) 


' That this formula affords a fair approximation even when 
m is small is seen from the case of a die thrown 12 times. The 
probability that any particular face will appear in one throw is 
p = 1/6, whence g = 5/6; and the most probable result in 12 
throws is that in which the particular face appears twice and 
fails to appear ten times. The probability of this result com- 
puted from (3) is 0.309, while the exact formula (1) gives 0.296. 


The probability that the event will occur a number of times 


* This expression is due to Laplace, Théorie Analytique des Probabilités, 
See also De Morgan’s Calculus, pp. 600-604, 


a a, 
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comprised between (4 — a) and (uw + a) in x trials is evidently 
expressed by the sum of the terms in (p +)" for which the 
exponent of f has the specified range of values. Calling this 
probability A, putting 


aa=nuptu, and m=ng—t, 


and using Stirling’s theorem (which implies that is a large 
number),* 


R 1 S I (1 mn uU ae : U lees 
~~ 2nnpq np ( ng ; 


very nearly; and the summation is with respect to w from 
u=—ato“z—=-+a. But expansion shows that the natural 
logarithm of the product of the two binomial factors in this 
equation is approximately — u’/2npg. Hence 


5 mipte 


= 
2 V20npq 


and, since 2 is supposed large, this may be replaced by a definite 
integral, putting 


dz—1/V2npqg, and 2 = u'/2npg. 


Thus 
Bay Vinpa a/ V2npq 
I 2 
a 0 a2 = e7*' d2. (4) 
V1 V1 if 
—a/V2npq M 


This equation expresses the theorem of James Bernoulli, 
given in his Ars Conjectandi, published in 1713. 


The value of the right-hand member of (4) varies, as it 
should, between o and 1, and approaches the latter limit rap- 
idly as 2 increases. Thus, writing for brevity 


Zz 


rag fords 


0 


* See Bertrand, Calcul des Probabilités, Paris, 1889, for an extended discus- 


’ sion of the questions considered in this Article. 
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the following table shows the march of the integral : 


z £ z F Zz og 
0.00 0.000 0.75 0.711 1.50 0.966 

25 276 I.00 843 1.75 987 

50 520 1.25 -923 2.00 995 


To illustrate the use of (4), suppose there is required the 
probability that in 6000 throws of a die the ace will appear a 
number of times which shall be greater than 1/6 % 6000 — IO 
and less than 1/6 X 6000-+ IO, or a number of times lying 
between ggo and 1o10. Inthis case a = 10, z = 6000, p = 1/6, 
qg=5/6. Thus, a/V2upg = 10/¥2.6000.1/6. 5/6 = 0.245. 
Hence, by (4) and the table, R = 0.27. 


Prob. 11. If the ratio of males to females at birth is 105 to 100, 
what is the probability that in the next 10,000 births the number of 
males will fall within two per cent of the most probable number ? 


Prob. 12. If the chance is even for head and tail in tossing a 
coin, what is the probability that in a million throws the difference 
between heads and tails will exceed 1500? 


ART. 7. INVERSE PROBABILITIES.* 


If an observed event can be attributed to any one of several 
causes, what is the probability that any particular one of these 
causes produced the event? To put the question in a concrete 
form, suppose a white ball has been drawn from one of two 
urns, U, containing 3 white and 5 black balls, and U, contain- 
ing 2 white and 4 black balls; and that the probability in favor 
of each urn is required. If U, is as likely to have been 
chosen as U,, the probability that U, was chosen is 1/2. After 
such choice the probability of drawing a white ball from UO is 
3/8. Before drawing, therefore, the probability of getting a 
white ball from U, was 1/2 X 3/8 = 3/16, by Art. 5. Similarly, 
before drawing the probability of getting a white ball from UC; 
was 1/2 X 2/6=1/6. These probabilities will remain un- 
changed if the number of balls in either urn be increased or 


* See Poisson, Probabilité des Jugements, pp. 81-83. 
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diminished so long as the ratio of white to black balls is kept 
constant. Make these numbers the same for the two urns. 
Thus let the first contain 9 white and 15 black, and the second 
8 white and 16 black; whence the above probabilities may be 
written 1/2 x 9/24 and 1/2 X 8/24. It is now seen that there 
are (9 + 8) cases favorable to the production of a white ball, 
each of which has the same antecedent probability, namely, 1/2. 
Since the fact that a white ball was drawn excludes considera- 
tion of the black balls, the probability that the white ball came 
from U, is 9/17 and that it came from U, is 8/17; and the sum 
of these is unity, as it should be. 

To generalize this result, let there be m causes, C,, C,,... Cy. 
Denote their direct probabilities by 7,, 9.,- ++ Qms their antecedent 
probabilities by 7,, 7,,--+7%mi and their resultant probabilities 
on the supposition of separate existence by Zio fas ss «La: 
That is, 

Pr=QUaty Pr= Vala e + Pm = Tul me (1) 

Let D be the common denominator of the right-hand mem- 
bers in (1), and denote the corresponding numerators of the 
several fractions by S,, 5,,.-+ 5. Then 

2,=5/D, p= Ue) aera pee ae ee Te 
and it is seen that there are in all (s;+5,+... Sm) equally 
possible cases, and that of these s, are favorable to C,, s, to 
ey peetience, 1 POL. s Fy denote the probabilities of 
the several causes on the supposition of their coexistence, 
P.=s f(y +5,4+ +++ $m) =P/(P, +2: 1+ + + Pun) 

Thus in general 

P.=7,/2p, Py = b,/ 2p, +++ Pu = bu/ 2P- (2) 

To illustrate the meaning of these formulas by the above 

concrete case of the urns it suffices to observe that 
for U., 9g,= 3/8 and 717,= 1/2, 
for U,, g,=1/3 and 7, = 1/2; 
whence p, = 3/16, pp = 1/6, py hr = 17/485 
and P.=90/17, P. = 8/17. 


Asa second illustration, suppose it is known that a white 
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ball has been drawn from an urn which originally contained m 
balls, some of them being black, if all are not white. What is 
the probability that the urn contained exactly white balls? 
The facts are consistent with s different and equally probable 
hypotheses (or causes), namely, that there were 1 white and 
(7 — 1) black balls, 2 white and (sz — 2) black balls, etc. 
Heneq:itt(1); op 9p = <= Sand 
PEE hy of = 2 iss Oh Se ie 

Thus =p = (1/2)(m-+ 1), 

2n 
He ee PL meme +- 1)° 

This shows, as it evidently should, that z = m is the most 
probable number of white balls in the urn. The probability 
for this number is P,, = 2/( + 1), which reduces, as it ought, 
toil. for 77 = 1. 

Formulas (1) and (2) may also be applied to the problem of 
estimating the probability of the occurrence of an event from 
the concurrent testimony of several witnesses, X,, X,,... 
Denote the probabilities that the witnesses tell the truth by 
4, ¥,,,.. Then, supposing them to testify independently, 
the probability that they will concur in the truth concerning 
the event is 7,4,...; while the probability that they will con- 
cur in the only other alternative, falsehood, is (I—#,)(I—+,)... 
The two alternatives are equally possible. Hence by equations 
(1) and (2) 


Pi = F440 65, Pee FA) = es, 


ns 

Fite ge (ie al Fee 

roe (he eee eee (3) 
ame freee ay Stee’ lee 3) bm 


P, being the probability for and P, that against the event. 


To illustrate (3), if the chances are 3 to 1 that X, tells the 
truth and 5 to 1 that X, tells the truth, x, = 3/4, #, = 5/O.and 
P,= 15/16; or, the chances are 15 to 1 that an event occurred 
if they agree in asserting that it did.* 


* For some interesting applications of equations (3) see note E of Appendix 
to the Ninth Bridgewater Treatise by Charles Babbage (London, 1838). 
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It is of theoretical interest to observe that if x,, 7,...in 
(3) are each greater than 1/2, P, approaches unity as the 
number of witnesses is indefinitely increased. 

Prob. 13. The groups of numbers of one figure each, two figures 
each, three figures each, etc., which it is possible to form from the 
nine digits 1, 2,...g9 are printed on cards and placed severally in 
nine similar urns. What is the probability that the number 777 will 
be drawn in a single trial by a person unaware of the contents of 
the urns? 


Prob. 14. How many witnesses whose credibilities are each 3/4 
are essential to make ?, = 0.999 in equation (3) ? 


ArT. 8. PROBABILITIES OF FUTURE EVENTS. 


Equations (2) of Art. 7 may be written in the following 
manner: 


ee (1) 

(iit ee epee 
If f,, Por + +» fm are found by observation, P,, /,,. . . P,, will ex- 
press the probabilities of the corresponding causes or their 
effects. When, as in the case of most physical facts, the num- 
ber of causes and events is indefinitely great, the value of any 
p or Pin (1) becomes indefinitely small, and the value of =f 
must be expressed by means of a definite integral. Let 4 de- 
note the probability of any particular cause, or of the event to 
which it gives rise. Then, supposing this and all the other 
causes mutually exclusive, (1 — x) will be the probability 
against the event. Now suppose it has been observed that in 
(m+ 2) cases the event in question has occurred 7 times and 
failed 2 times. The probability of such a concurrence is, by 
Art. 5,cx(1 — x)", where c is a constant. Since x is unknown, 
it may be assumed to have any value within the limits o and 1; 
and all such values are a priori equally possible. Put 

p= te"0.— 4)". 

Then evidently the probability that + will fall within any as- 
signed possible limits @ and @ is expressed by the fraction 


[vex] fvde: 
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so that the probability of any particular + is given by 


a1 — 4)"dx 
ps "(1 — 4r)"dx 


This may be regarded as the antecedent probability of the 
cause or event in question. 


P= (2) 


What then is the probability that in the next (7 + s) trials 
the event will occur 7 times and fail s times, if no regard is had 
of the order of occurrence? If x were known, the answer 
would be by Arts, 2 and 5 

(aes) ae , 

Bee: Bersted (3) 

But since 4 is restricted only by the condition (2), the required 

probability will be found by taking the product of (2) and (3) 

and integrating throughout the range of x. Thus, calling the 
required probability Q, 


1 
Kt = 2y¥a¢ 


O = (r+s)! o 
: ris! } : (4) 
JES — s)"dx 
The definite integrals which appear here are known as Gamma 
functions. They are discussed in all of the higher treatises on 
the Integral Calculus. Applying the rules derived in such 
treatises there results * 


_&r+s)!\(m@+r'(n+5)!\(m+n+1)! 
— ris!minl(m+nutr+ts +1)! * (5) 
If regard is had to the order of occurrence of the event ; 


that is, if the probability required is that of the event happen- 
ing 7 times in succession and then failing s times in succession, 


* It isa remarkable fact that formula (5) is true without restriction as to 
values of m, »,r, s. The formula may be established by elementary considera- 
tions, as was done by Prevost and Lhuilier, 1795. See Todhunter’s History of 
he Theory of Probability, pp. 453-457. 
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the factor (r+ s)!/r!s! in (3), (4), (5) must be replaced by 
unity. 

To illustrate these formulas, suppose first that the event 
has happened mm times and failed no times. What is the prob- 
ability that it will occur at the next trial? In this case (4) 
gives 


= furrax | fvrax = (m+ 1)/(m 4-2). 


When wm is large this probability is nearly unity. Thus, the 
sun has risen without failure a great number of times m; the 
probability that it will rise to-morrow is 


Cope eee 


m m ne eee 
which is practically 1. 


Secondly, suppose an urn contains white and black balls in 
an unknown ratio. If in ten trials 7 white and 3 black balls 
are drawn, what is the probability that in the next five trials 
2 white and 3 black balls will be drawn? The application of 
(5) supposes the ratio of the white and black balls in the urn 
to remain constant. This will follow if the balls are replaced 
after each drawing, or if the number of balls in the urn is sup- 
posed infinite. The data give 


Bie ae F025 $= 3, 
mtr=9 u+s=6 rvt+s=5, mt+n+1=I11, 
mtutrts+i1= 16. 
Thus by (5) 
s!o!6!rr! 


Olas igliela (ok 


Suppose there are two mutually exclusive events, the first 
of which has happened m times and the second z times in 
m-+ntrials. What is the probability that the chance of the 
occurrence of the first exceeds 1/2? The answer to this ques- 
tion is given directly by equation (2) by integrating the nume- 
rator between the specified limits of x. That is, 


a 
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] 


iy a1 — £)"ax 


5 


eee (6) 
i, ‘e"(1 — 4«)"dx 


Thus, if #—= 1 and % =o, P= 3/4; or the odds are three to 
one that the event is more likely to happen than not. Simi- 
larly, if the event has occurred: 7 times in succession, 


Pal A(t 2h, 


which approaches unity rapidly with increase of 2. 


ART. 9. THEORY OF ERRORS. 


The theory of errors may be defined as that branch of math- 
ematics which is concerned, first, with the expression of the re- 
sultant effect of one or more sources of error to which com- 
puted and observed quantities are subject ; and, secondly, with 
the determination of the relation between the magnitude of 
an error and the probability of its occurrence. In the case of 
computed quantities which depend on numerical data, such as 
tables of logarithms, trigonometric functions, etc., it is usually 
possible to ascertain the actual values of the resultant errors. 
In the case of observed quantities, on the other hand, it is not 
generally possible to evaluate the resultant actual error, since 
the actual errors of observation are usually unknown. In either 
case, however, it is always possible to write down a symbolical 
expression which will show how different sources of error enter 
and affect the aggregate error; and the statement of such an 
expression is of fundamental importance in the theory of errors. 


To fix the ideas, suppose a quantity Q to be a function of 
several independent quantities x, y, z...; that is, 


O= (ey reare), 
and let it be required to determine the error in Q due to errors 
in 4,-¥, 2... Denote such errors by 40,427, Aye 
Then, supposing the errors so small that their squares, prod- 
ucts, and higher powers may be neglected, Taylor’s series gives 
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AO i ar teh das a 48s: (1) 


This equation may be said to express the resultant actual error 
of the function in terms of the component actual errors, since 
the actual value of 4Q is known when the actual errors of 
4,7, 2...are known. It should be carefully noted that the 
quantities 4, y, Z... are supposed subject to errors which are 
independent of one another. The discovery of the independent 
sources of error is sometimes a matter of difficulty, and in general 
requires close attention on the part of the student if he would 
avoid blunders and misconceptions. Every investigator in work 
of precision should have a clear notion of the error-equation of 
the type (1) appertaining to his work; for it is thus only that 
he can distinguish between the important and unimportant 


sources of error. 


Prob. 15. Write out the error-equation in accordance with (r) 
for the function Q = xyz + * log (y/2). 

Prob. 16. In a plane triangle a/d = sin A/sin B. Find the error 
in a due to errors in 4, 4, and B. 

Prob. 17. Suppose in place of the data of problem 16 that the 
angles used in computation are given by the following equations : 
A=4,+4(180°—A,— B,— C,), B= B,+4(180° — 4,—8, — C,), 
where 4,, B,, C, are observed values. What then is 4a? 

Prob. 18. If w denote the weight of a body andr the radius of 
the earth, show that for small changes in altitude, 4z/w=—24r/r; 
whence, if a precision of one part in 509 000 000 is attainable in com- 
paring two nearly equal masses, the effect of a difference in altitude 
of one centimeter in the scale-pans of a balance will be noticeable.* 


Art. 10. Laws OF ERROR. 


A law of error is a function which expresses the relative 
frequency of occurrence of errors in terms of their magnitudes. 
Thus, using the customary notation, let e denote the magni- 


* This problem arose with the International Bureau of Weights and Measures, 
whose work of intercomparison of the Prototype Kilogrammes attained a pre- 
cision indicated by a probable error of 1/500 000 oooth part of a kilogramme. 
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tude of any error ina system of possible errors. Then the law 
of such system may be expressed by an equation of the form 
I = $(€)- (1) 
Representing € as abscissa and y as ordinate, this equation 
gives acurve called the curve of frequency, the nature of which, 
as is evident, depends on the form of the function @ This 
equation gives the relative frequency of occurrence of errors in 
the system ; so that if € is continuous the probability of the 
occurrence of any particular error is expressed by yde = (€)de; 
which is infinitesimal, as it plainly should be, since in any con- 
tinuous system the number of different values of € is infinite, 
Consider the simplest form of ¢(e), namely, that in which 
@(€) =c¢,aconstant. This form of d(€) obtains in the case of 
the errors of tabular logarithms, natural trigonometric func- 
tions, etc. In this case all errors between minus a half-unit 
and plus a half-unit of the last tabular place are equally likely 
to occur. Suppose, to cover the class of cases to which that 
just cited belongs, all errors between the limits —@and +a 
are equally likely to occur. The probability of any individual 
error will then be ¢(e)de = cde, and the sum of all such prob- 
abilities, by equation (2), Art. 4, must be unity. That is, 


+a +a 
Jf elede =e fide =1. (2) 
This gives ¢ = 1/2a, or by (1) y= 1/2a. The curve of fre- 
Q quency in this case is shown in the figure, 
C) p AB being the axis of € and OQ that of y. 


It is evident from this diagram that if the 

errors of the system be considered with 

respect to magnitude only, half of them 

should be greater and half less than @/2. 

a 5 se This is easily found to be so in the case of 
tabular logarithms, ete. 

As a second illustration of (1), suppose y and € connected 


by the relation y = ¢ Va’ — e’, where a is the radius of acircle, 
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c a constant, and € may have any value between — aand +a. 


Then the condition 
+a 


c yi deVa—€=1 
gives c = 2/(a’a). In this, as in the preceding case, (+ €) = 
¢(— e), the meaning of which is that positive and negative 
errors of the same magnitude are equally likely to occur. It 
will be noticed, however, that in the latter case small errors 
have a much higher probability than those near the limit a, 
while in the former case all errors have the same probability. 


In general, when € is continuous ¢(¢e) must satisfy the condi- 
tion if, d(e)de = 1, the limits being such as to cover the entire 


range of values of «. The cases most commonly met with are 
those in which @(é) is an even function, or those in which 
d(+6=¢(— 4. In such cases, if + @ denote the limiting 
value of «€, 


[dene = 2 f ede ats (3) 


Art. 11. TYPICAL ERRORS OF A SYSTEM. 


Certain typical errors of a system have received special 
designations and are of constant use in the literature of the 
theory of errors. These special errors are the probable error, 
the mean error, and the average error. The first is that error 
of the system of errors which is as likely to be exceeded as 
not ; the second is the square root of the mean of the squares 
of all the errors; and the third is the mean of all the errors 
regardless of their signs. Confining attention to systems in 
which positive and negative errors of the same magnitude 
are equally probable, these typical errors are defined mathe- 
matically as follows. Let 

é, = the probable error, 
€,== the mean error, 
"e, = the average error. 
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Then, observing (2) of Art. 10, 


—e +e 
df werde = =/ pe)de = wede = =) pede = 4 | 


in, ie) 
En = | pejerde, = 2 fo eyde | 


The student should seek to avoid the very common misap- 
prehension of the meaning Ae the probable error. It is not 
“the most probable error,” nor “the most probable value of 
the actual error’; but it is nee error which, disregarding signs, 
would occupy the middle place if all the errors of the system 
were arranged in order of magnitude. A few illustrations will 
suffice to fix the ideas as to the typical errors. Thus, take the 
simple case wherein #(€) =¢ = 1/2a, which applies to tabular 
logarithms, etc. Equations (1) give at once 


I == 
= tb 54 én = ONS, é= oa. 


For the case of tabular values, a= 0.5 in units of the last 
tabular place. Hence for such values 
Gy = 2b. 0.28, y= 0.20, pene 105. 

Prob. 19. Find the typical errors for the cases in which the law 
of error is P(e) = cVa? — &, fle) =e 4aFe), P(e) =c cos*(ze/2a); 
c being a constant to be determined in each case and € having any 
value between — a and + @. 


ART. 12. LAWS OF RESULTANT ERROR. 


When several independent sources of error conspire to pro- 
duce a resultant error, as specified by equation (1) of Art. 9, 
there is presented the problem of determining the law of the 
resultant error by means of the laws of the component errors. 
The algebraic statement of this problem is obtained as follows 
for the case of ‘continuous errors: 

In the equation (1), Art. 9, write for brevity 


md _ ad 
=a, €é,= pe é,=— —d4y,. oe 


LAWS OF RESULTANT ERROR, SY) 


and let the laws of error of €, €,, €,,...- be denoted by ¢(e), 
p,(e,), ,(€,).-. Then the value of e is given by 
e=e,te,+... (1) 
The probabilities of the occurrence of any particular values 
of €, €,,...are given by @(€,)de,, P,(&)e, - + - 5 and the 
probability of their concurrence is the probability of the cor- 
responding value of ¢. But since this value may arise in an 
infinite number of ways through the variations of €,, €,,... 
over their ranges, the probability of e, or $(e)de, will be 
expressed by the integral of p,(e,)de,p,(€,)@e, . .. subject to 
the restriction (1). This latter gives e, = e— € — +++, and 
de, = de for the multiple integration with respect to ¢,, €,,< . 
Hence there results 


d(ej)de = de [ H(¢ —e,— €,—...)P,(€,)de,-- +, 


or 
P(€) = fole-4 Sy SO Jole,)de, f o.(es)4es cust? (2) 


It is readily seen that this formula will increase rapidly in 
complexity with the number of independent sources Disctror.= 
For somé of the most important practical applications, how- 
ever, it suffices to limit equation (2) to the case of two inde- 
pendent sources of error, each of constant probability within 
assigned limits. Thus, to consider this case, let €, vary over 
the range — a to +a, and e¢, vary over the range — 6to+ 0. 
Then by equation (2), Art. 10, 

o,(€,) = 1/(2a),  ,(€,) = 1/(24). 


Hence equation (2) becomes 


I 
f(€) = zai, 1 


In evaluating this integral «, must not surpass +6 and 
€, = € — €, must not surpass + @. Assuming a > 4, the limits 
of the integral for any value of €= eé, te, lying between 
— (a+ 6) and —(@ — 4) are — band +(e-+.a). This fact is 

* The reader desirous of pursuing this phase of the subject should consult 


Bessel’s Untersuchungen ueber die Wahrscheinlichkeit der Beobachtungsfehler; 
Abhandlungen von Bessel (Leipzig, 1876), Vol. II. 
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made plain by a numerical example. For instance, suppose 
a=s5andé=3. Then — (a+ 6) = —8and — (a— 4) = — 3. 
Take e = — 6,a number intermediate to — 8 and — 3. Then 
the following are the possible integer values of e, and €, which 
will produce e = — 6: 


€ &; are, limits of e, 
—6=—5-—1, —1=+(e+ 4), 
ea A as 


=> d 
Similarly, the limits of e, for values of e lying between 
— (a — 6) and + (a@— 8) are — 6 and + 4; and the limits of 
e, for values of e between + (a—6) and + (a+ 4) are +(€—a) 
and +06, Hence 


e+a 
Ho=si5 | t= cae for —(a+d) <e< — (a2—4), 


—b 


HO= Fy f dyes for end) ce  aeee G3) 
40 
6)= TF fe eae for +(a—b)<e<+(a+0). 


€-a “J 


Thus it appears that in this case the graph of the resultant 
law of error is represented by the upper base and the two sides 
of a trapezoid, the lower base being the 


IIOII ; ; re ; 

AE axis of e and the line joining the middle 
oint s bei : 

Henin points of the bases being the axis of P(E). 


(See the first figure in Art. 13.) The prop- 
erties of (3), including the determination 
of the limits, are also illustrated by the 
adjacent trapezoid of numerals arranged 
to represent the case wherein a = 0.5 and 
6=0.3. The vertical scale, or that for P(e), does not, how- 
ever, conform exactly to that for e. 


TTL EOI iit 
Pe TOU wes 
IITIIIIIOIVIIIIII 
IITILIIIIOTIIIIIII 
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Prob. 20. Prove that the values of #(€) as given by equation (3) 
satisfy the condition specified in equation (3), Art. 10. 

Prob. 21. Examine equations (3) for the cases wherein a = é and 
b = o; and interpret for the latter case the first and last of (3). 

Prob. 22. Find from (3), and (1) of Art. 11, the probable error of 
the sum of two tabular logarithms. 


Art. 13. ERRORS OF INTERPOLATED VALUES. 


Case I.—One of the most instructive cases to which formulas 
(3) of Art. 12 are applicable is that of interpolated logarithms, 
trigonometric functions, etc., dependent on first differences. 
Thus, suppose that v, and v, are two tabular logarithms, and 
that it is required to get a value v lying 7 tenths of the interval 
from v, towards v, Evidently 

v=r,ty, —v,)f=(U — 2), + %5 
and hence if ¢, @, ¢, denote the actual errors of ¥, v,, , , fe- 
spectively, 
e=(1 —2e,-+ %e,. (1) 

It is to be carefully noted here that ¢as given by (1) re- 
quires the retention in v of at least one decimal place be- 
yond the last tabular place. For example, let v = log (24373) 
from a 5-place. table. Then %,= 4.38686, v, = 4.38703, 
v7, — V, = + 0.00017, = 0.3, and v = 4.38691.1. Likewise, as 
found from a 7-place table, ¢, = — 0.45, @& = + 0.37 in units of 
the fifth place; and hence by (1)¢= —0.20. That is, the 
actual error of v = 4.38691.1 is = 0.20, and this is verified by 
reference to a 7-place table. 

The reader is also cautioned against mistaking the species 
of interpolated values here considered for the species common- 
ly used by computers, namely, that in which the interpolated 
value is rounded to the nearest unit of the last tabular place. 
The latter species is discussed under Case II below. 

- Confining attention now to the class of errors specified by 
equation (1), there result in the notation of the preceding 


article 
6=(1—Z)e, & = le and ¢=¢ =, + €,; 


and since e, and e, each vary continuously between the limits 
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+0.5 of a unit of the last tabular place, a and 6 in equations 
(3) of that article have the values 
a= 0.5(1 = 2). 0 — Gabe 
Hence the law of error of the interpolated values is ex- 
pressed as follows: 


9(€).= a —_ values of € betw. —o.5 and —(0.5 —2), 


it as 
ON — & 


— ———- for values of € betw. +(0.5—2) and +0.5. 
(=o 


| 
= =~ for values of ' betw. —(0.5—¢) and +(0.5 —2), | (2) 
J 
The graph of #(€) for ¢=1/3 is shown by the trapezoid 
AB, BC, CD in the figure on page 40. Evidently the equa- 
tions (2) are in general represented by a trapezoid, which degen- 
erates to an isosceles triangle when ¢ = 1/2. 
The probable, mean, and average errors of an interpolated 
value of the kind in question are readily found from (2), and 
from equations (1) of Art. 11, to be 


e = (/4ti 2) for oO faa a) 
= 1/2 — (1/2)Vari— 7) for 1/37 2 2/3 
= s/t fori2/ 3) <7 mal; 

1 — (1 — 2)? | 
HES) 96(1 — Az r (3) 
6. = oe for) 0 <f aan 2, 

ml ST fort (2h. | 


It is thus seen that the probable error of the interpolated 
value here considered decreases from 6.25 to 0.15 of a unit of 
the last tabular place as ¢ increases from o to 0.5. Hence such 
values are more precise than tabular values; and the computer 
who desires to secure the highest attainable precision with a 
given table of logarithms should retain one additional figure 
beyond the last tabular place in interpolated values. 
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Case II.—Recurring to the equation v = v, + “(v, — v,) for an 
interpolated value v in terms of two consecutive tabular values 
v, and v,, it will be observed that if the quantity tv, — v,) 1s 
rounded to the nearest unit of the last tabular place, a new error 
is introduced. For example, if v,=log 1633 = 3.21299, and 
v, = log 1634 = 3.21325 from a 5-place table, v, — v,= + 26 
units of the last tabular place; and if ¢= 1/3, ¢(v, — v,) = 83; 
so that by the method of interpolation in question there results 
v = 3.21299 + 9 = 3.21308. Now the actual errors of v, and 
vy, are, as found froma 7-place table, — 0.38 and + 0 21 in units 
of the fifth place. Hence the actual error of v is by equation 
(1), #X — 0.38-+ 3 X + 0.21 — 4 = — 0.52, as is shown di- 
rectly by a 7-place table. 

It appears, then, that in this case the error-equation cor- 
responding to (1) is 

ea(r— ae, +t, +4, ee 
wherein e, and e, are the same as in (1) and e, is the actual'error 
that comes from rounding ¢(v, — v,) to the nearest unit of the 
last tabular place. 

The error ¢,, however, differs radically in kind from e, and 
é, The two latter are continuous, that is, they may each have 
any value, between the limits — 0.5 and + 0.5; while e, is dis- 
continuous, being limited to a finite number of values depend- 
ent on the interpolating factor 4 Thus, for f= 1/2 the only 
possible values of ¢, are 0 4. 1/2, and — 1/2; likewise for ¢ = 
1/3, the only possible values of e, are 0, + 1/3, and — 1/3. It 
is also clear that the maximum value of ¢, which is constant and 
equal to 1/2 for (1), is variable for (4) in a manner dependent 
onz. For example, in (4); 

The maximum of e=1/2-+1/2=1, fort= 1/2) 
7 ‘. “e@=1/2 +1/3= 5/6,“ t= 1/3) 
a ce g—1/2+1/2=1,° “ =i A, 
a am “e=1/2+ 2/5 =9/10" ft 5. 

The determination of the law of error for this case presents 
some novelty, since it is essential to combine the continuous 
meirors (1 — fe, and Ze, with the discontinuous error ¢, The 
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simplest mode of attacking the problem seems to be the fol- 
lowing quasi-geometrical one. In the notation of Arts. 12 and 
13, put in (4) e=e, (1—Z)e, = €,, fe, =e, and¢,=e,. Then 

e = (€,+ €,)-+ & (5) 

The law of error for (€, + €,) is given by equation (2) for any 
value of ¢. Hence for a given value of ¢ there will be as many 
expressions of @(e) as there are different values of €,. The 
graphs of (e) will all be of the same form but will be differently 
placed with reference to the axis of P(e). Thus, if ¢= 1/3 the 

values of e, are — 1/3, O, and 

+ 1/3, and these are equally 

likely tO dccur. FOr ep) ¢ Ale 

graph is given directly by (2), 

and is the trapezoid ABCD 

symmetrical with respect to OQ. 
For ¢,=—1/3 thé graph is 
\ abQd, of the same form as 
\' ABCD but shifted to the left 
by the amount of e,= — 1/3. 
aq b a Similarly, the graph for the case 
of e, = + 1/3 is a’Qd’a’, and is produced by shifting ABCD to 
the right by an amount equal to + 1/3. 

Now, since the three systems of errors for this case are 
equally likely to occur, they may be combined into one system 
by simple addition of the corresponding element areas of the 
several graphs. Inspection of the diagram shows* that the 
_resultant law of error is expressed by 


He) = (1/4)(5 + 6e) for —s/6<e<—1/6, 
=1 for —1/6 <e<+ 1/6, (6) 
= (1/4)(5 —6e) for+1/6<e<c4+5/6, 


This is represented by a trapezoid whose lower base is 10/6, 
upper base 2/6, and altitude 1. 


* Sum the three areas and divide by 3 to make resultant area — I 
4 , 


as 
required by equation (3), Art. 10. 
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As a second illustration, consider equation (5) for the case 
t= 1/2. In this case e, must be either 0 or 1/2, the sign of 
which latter is arbitrary. For e,= 0, equations (2) give 

ge) =2+4e for —1/2 <€ <0, 
=2-—4e for ae a (7) 
This function is represented by the isosceles triangle AQ 
whose altitude OQ is twice the base AZ. 

Similarly ¢(¢) for e, = + 1/2 would 
be represented by the triangle AQE dis 
placed to the right a distance 1/23 and 
if the two systems for ¢,=o0 and €, = 
+ 1/2 be combined into one system, 
their resultant law of error is evidently 
@(e)= 1+2€ for —1/2<€<0, 

== 1 for 0< €<--1/2, (8) 

= 9 %e for+-1/z<.€< 1; 4 
the graph of which is ABCD. On the 
other hand, if the errors in this combined system be considered 


O BD) D 


with respect to magnitude only, the law of error is 
ple 2(1-— ¢) tors O-<e<st, (9) 
the graph of which is OOD. 
The student should observe that (6), (7), (8), and (9) satisfy 


the condition J eed — 1 if the integration embraces the 


whole range of «. 

The determination of the general form of P(e) in terms of 
the interpolating factor ¢ for the present case presents some 
difficulties, and there does not appear to be any published solu- 
tion of this problem.* The results arising from one phase of 
the problem have been given, however, by the author in the 
Annals of Mathematics,t and may be here stated without 
proof. The phase in question is that wherein ¢ is of the form 
1/n, m being any positive integer less than twice the greatest 


® The author explained a general method of solution in a paper read at the 
summer meeting of the American Mathematical Society, August, 1895. 


¢ Vol. II, pp. 54-59. 
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tabular difference of the table to which the formulas are ap- 
plied. For this restricted form of ¢ the possible maximum 
value of € as given by equation (5) is, in units of the last 
tabular place, (2% — 1)/n for m odd and 1 for even, 

The possible values of €, of equation (5) are 


I 2 i — 1 


oO — — for z odd 
» ot n’ L n? 2n : 
a I a 2 ” 2D i I PF 
Oo = —,...-+- — for 2 even. 
: n°’ n’ ON 2 


An important fact with regard to the error 1/2 for z even 
is that its sign is arbitrary, or is not fixed by the computation 
as is the case with all the other errors. However, the com- 
puter’s rule, which makes the rounded last figure of an inter- 
polated value even when half a unit is to be disposed of, will, 


in the long-run, make this error as often plus as minus. 


The laws of error which result are then as follows: 


For 2 odd. 
p(e) = 1 for e between — 1/2z and + 1/22, 
P(e) = 4 feet +P e) for ebetw. -1/2z and F(2n—1)/2n. 
na—I1\ 2n 
For 2 even. 
n ZIU—— 22 
P(e) = Go a + e) for e between 0 and F 1/n, 
nN Zita 
= a a + e) for € betw. F 1/z and F (wm — 1)/n, 
= Fae eae + €) for e between (n— 1/m) and Ft. 


By means of these formulas and (1) of Art. 11 the probable, 
mean, and average errors for any value of can be readily 
found. The following table contains the results of sucha com- 
putation for values of z ranging from 1 to 10. The maximum 
actual error for each value of z is also added. The verifica- 


tion of the tabular quantities will afford a useful exercise to the 
student. 


ERRORS OF INTERPOLATED VALUES. 43 


TypicAL Errors OF INTERPOLATED LOGARITHMS, ETC. 


Interpotating Probable Error, Mean Error. | Average Error, : 
Factor. € € € Maximum 
t=1/n b m a Actual Error. 
I 0.250 0.289 0.250 1/2 
1/2 .292 .408 +333 I 
1/3 +256 347 .287 5/6 
1/4 .276 382 +313 I 
1/5 -208 +370 303 9/10 
a 277 ae +315 I 
A +274 +300 “GEL 13/14 
1/3 .279 .389 318 I 
1/9 .278 -386 .316 17/18 
1/10 .281 -392 -320 I 


When the interpolating factor ¢ has the more general form 
ne /It, wherein m and are integers with no common factor, the 
possible values of €, are the same as forf=—1/a. But equa- 
tions (3) of Art. 12 are not the same for ¢= m/n as fort = 1/n, 
and hence for the more ceneral form of 7, @(e€) assumes a new 
type which is somewhat more complex than that discussed 
above. The limits of this work render it impossible to extend 
the investigation to these more complex forms of #(€). It may 
suffice, therefore, to give a single instance of such a function, 
namely, that for which ¢ = 2/5. For this case 

o(6) = 1 for e between 0 and F 1/10, 

(5 /6)(13/10 + e) fore between — 1/10 and F 3/10, 
(5/3)(4/5 + €) fore between # 3/10 and =F 7/10, 
— (5/6)(9/1o+ €) for e between + 7/10 and = 9Q/IO. 

The graph of the right-hand half of 4 8 
this function is shown in the accompany- 
ing diagram, the whole graph being 
symmetrical with respect to OA, or the 
axis of P(€). 

_ Attention may be called to the strik- 
ing resemblance of this graph to that of 
the law of error of least squares. 0 


lI 


I 


Prob. 23. Show from equations (3) that ¢,, varies from 1/V1 
= 0.29 —, for = 09, to 1/V¥ 24 =0.20-+, for = 0.5; and that €, 
varies from 0.25 to 1/6 for the same limits. 
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Prob. 24. Show that the probable, mean, and average errors 
for the case of ¢= 2/5 cited above (p. 43) are + 0.261, + 0.251, 
and + 0.290, respectively. 


ART. 14. STATISTICAL TEST OF THEORY. 

A statistical test of the theory developed in Art. 13 may 
be readily drawn from any considerable number of actual er- 
rors of interpolated values dependent on the same interpolating 
factor. The application of such a test, if carried out fully by 
the student, will go far also towards fixing clear notions as to 
the meaning of the critical errors. 

Consider first the case in which an interpolated value falls 
midway between two consecutive values, and suppose this 
interpolated value retains two additional figures beyond the 
last tabular place. Then by equations (2), Art. 13, the law of 
error of this interpolated value is 

o(e) = 2-+ 4e for e between — 0.5 ando 
= 2 — 4e for e between 0 and+- 0.5. 


Hence by equation (1) of Art. 11, or equation (3) of Art. 12, the 
probable error in this system of errors is $ — (1)¥2=0.15. 
It follows, therefore, that in any large number of actual errors 
of this system, half should be less and half greater than 0.15. 
Similarly, of the whole number of such errors the percentage 


falling between the values 0.0 and 0.2 should be 
+o0-2 +o-2 


JS $Qde = = af (2 — 4€)de = 0.64; 


-0.2 
that is, sixty-four per cent of the errors in question should be 
less numerically than 0.2. 

To afford a more detailed comparison in this case, the act- 
ual errors of five hundred interpolated values from a 5-place 
table have been computed by means of a 7-place table. The 
arguments used were the following numbers: 20005, 20035, 
20065, 20105, 20135, etc., in the same order to 36635. The 
actual and theoretical percentages of the whole number of 
errors falling between the limits 0.0 and 0.1, 0.1 and 0.2, ete, 
are shown in the tabular form following: » ~ 
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oe rors Pare Pecconase 
Coram Ons st. ss vp Chie ee Verere nae ees 23:2 30 
Cleat One yaa ce ces sites s ots Pe eae) 28 
eta Os? firs s. < ie Swiss eile sims Zee 19.0 - 20 
Cn Osx oe ec Sik sss weale, seh CALS om We) I2 
ACAI Oss finns ses < ss SOS SURG AO Se 4.4. 4 
OEO. Dy. cass cavasescs bas on coke led 50 


The agreement shown here between the actual and theoretical 
percentages is quite close, the maximum discrepancy being 2.8 
and the average 1.5 per cent. 

Secondly, consider the case of interpolated mid-values of the 
species treated under Case II of Art. 13. The law of error for 
this case is given by the single equation (9) of Art. 13, namely, 
P(e) = 2(1 — e), no regard being paid to the signs of the errors. 
The probable error is then found from 


&p 
2 f(r — ede = , 


whence €, = 1 —$¥V2=0.29. Similarly, the percentage of 
the whole number of errors which may be expected to lie, for 
example, between 0.0 and ©.2 in this system is 


2 f(s — e)de = 0.36. 


Using the same five hundred interpolated values cited 
above, but rounding them to the nearest unit of the last tabu- 
lar place and computing their actual errors by means of a 7-place 
table, the following comparison is afforded : 


A Actual Theoretical 
sa ace ao te Percentage, Percentage. 
OOO atee so wo ole oiera son eee wheres 35.0 36 
Ris ACO ae he ase oe ipenieeies meee 2700 28 
OP ANCLOO tnd sud a's irs pp eg gota eed 150 20 
OOO Ob.cae asa <r henaavinn nace Lese 12 
SP Os wale sak ease ease? $0 “4 


O10 aNd 0.29.0 66 Sac code smite Ty -4Q.8 50 
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The agreement shown here between the actual and theoretical 
percentages is somewhat closer than in the preceding case, the 
maximum discrepancy being only 1.6 and the average only 0.6 
per cent. 

Finally, the following data derived from one thousand act- 
ual errors may be cited. The errors of one hundred inter- 
polated values rounded to the nearest unit of the last tabular 
place were computed * for each of the interpolating factors 
0.1, 0.2,..-0.9. The averages of these several groups of act- 
ual errors are given along with the corresponding theoretical 
errors in the parallel columns below: 


Interpolating Actual Theoretical 
Factor. Average Error. Average Error. 
OS Pig oe i OR ET TO LE 0.320 
Oana: Es BOLTS Peet POLS OL 0.288 0.303 
SE in OE eee IR Ci 0.321 0.304 
OiAte ane veils! «os Gh attal® Sick es IO: 208 0.290 
oS ee ceeats tas «6 see Sener ewer O324 0.333 
GASES eg OO a er Tk 0.276 0.290 
Oss Es Be Ak co ao ete O21 0.304 
O82 lee ota Fo ears Pied dete eews 0.289 0.303 
O.Ois 0 vised oo 8 Watbaetenaine ees bane Osa7 0.320 


The average discrepancy between the actual and theoret- 
ical values shown here is 0.017. It. is, perhaps, somewhat 
smaller than should be expected, since the computation of the 
actual errors to three places of decimals is hardly warranted 
by the assumption of dependence on first differences only. 


The average of the whole number of actual errors in this 
case is 0.308, which agrees to the same number of decimals 
with the average of the theoretical errors. + 


* By Prof. H. A. Howe. See Annals of Mathematics, Vol. III, p. 74. 
The theoretical averages were furnished to Prof. Howe by the author. 

+ The reader who is acquainted with the elements of the method of least 
squares will find it instructive to apply that method to equation (1), Art. 13, 
and derive the probable error of ¢. This is frequently done without reserve by 
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Prob. 25. Apply formulas (3) of Art. 12 to the case of the sum 
or difference of two tabular logarithms and derive the correspond- 
ing values of the probable, mean, and average errors. ‘The graph 
of £(€) is in this case an isosceles triangle whose base, or axis of €, 
is 2, and whose altitude, or axis of £(€), is 1. 


those familiar with least squares. Thus, the probable error of ¢ or e2 being 
0.25, the probable error of ¢ is found to be 


0.25 Vi — 27 + 27?, 


This varies between 0.25 for ¢=0 and 0.18 for = 3 while the true value of 
the probable error, as shown by equations (3), Art. 13, varies from 0.25 to 0.15 
for the same values of ¢. It is, indeed, remarkable that the method of least 
squares, which admits infinite values for the actual errors é1 and é, should give 
so close an approximate formula as the above for the probable error of e. 

Similarly, one accustomed to the method of least squares would be inclined 
to apply it to equation (4), Art. 13, to determine the probable error of e. The 
natural blunder in this case is to consider é1, é2, and és independent, and e; like 
1 and ¢, continuous betweer the limits 0.0 and 0.5; and to assign a probable 
error of 0.25 to each. In t’.is manner the value 


0.25 V2(1. — ¢-+ 2?) 


is derived. But this is absurd, since it gives 0.25 V2 instead of 0.25 forf=o0, 
The formula fails then to give even approximate results except for values of ¢ 
near 0.5, 
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